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^ , Abstract. The moduli space of principally polarized abelian surfaces is parametrized 

by three Igusa functions. In this article we investigate a new way to evaluate these 
functions by using Siegel Eisenstein series. We explain how to compute the Fourier 
coefficients of certain Siegel modular forms using classical modular forms of half- 
integral weight. One of the results in this paper is an explicit algorithm to evaluate 
the Igusa functions to a prescribed precision. 
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1. Introduction 



The classical theory of complex multiplication gives an explicit description of the 
Hilbert class field of an imaginary quadratic field: for a fundamental discriminant 
D < 0, the Hilbert class field of K = Q(>/D) is obtained by adjoining the value 
j((D + \AD)/2) to K. Here, j : H — >■ C is the classical modular function with 
Fourier expansion j(z) = 1/q + 744 + 196884^ + . . . in q = ex.p(2Triz). There are 
various ways to compute the minimal polynomial of j((D + vU)/2), and one of 
the most frequently used approaches proceeds by evaluating the j-function to high 
precision. 

The j-function is invariant under the action of SL2(Z) on the upper half plane H. 
To evaluate j(r), we may assume that r is in the 'standard' fundamental domain 
for SL2(Z)\H as described in e.g. [21, Sec. VILLI]. The naive approach to evaluate 
j(r) is to simply compute enough Fourier coefficients using for instance the recursive 
formulas given in [19]. Alternatively, one can use the relation 



j(z) = 1728 . .f ^j- . .- (1.1) 

g 2 (z) 3 - 27 g 3 (z) 2 

expressing the j-function in terms of the normalized Eisenstein series #2>#3 of 
weight 4 and 6. Better results can be obtained [1] by using the Dedekind 77-function 
defined by n(z) = a 1 ' 24 " n^Li(l — Q. n )i an d which satisfies 

_( V (z/2)/ V (z)r + W 
3\ z 



(rj(z/2)/rj(z)r 
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The sparsity of the g-expansion of the ry-function makes it very efficient for explicit 
computations. 

The j-function is intrinsically linked to the theory of elliptic curves, and the 
situation outlined above can be viewed as the '1-dimensional' case of complex mul- 
tiplication theory. In dimension 2, suitably chosen invariants of principally polarized 
abelian surfaces generate abelian extensions of degree 4 CM-fields, see [22, Sec. 15] 
for a precise statement. A popular choice of invariants are the three Igusa functions 
ji,J2,J3 defined below. Just as evaluating the elliptic j-function has applications 
to elliptic curve cryptography, evaluating Igusa functions is an important step in 
construction genus 2 curves suitable for use in cryptography, see e.g. [25]. 

The explicit evaluation of Igusa functions is less developed than its dimension- 1 
counterpart. Most people use ^-functions to evaluate Igusa functions. The (rather 
unwieldy) formulas expressing Igusa functions in terms of ^-functions are given in 
e.g. [25, pp. 441-442]. There is also a direct analogue of formula (1.1) which ex- 
presses the Igusa functions as rational functions in the Siegel Eisenstein series E w . 
Indeed, Igusa [11, p. 195] defines the normalized cusp forms 

43867 

Xio - ~2 12 -3 5 -5 2 -7-53 ( 4 6 ~ 10 > 



and 



* 12 = 2^ 31 53 59 72.337 (32 ' 72El + 2 ' ^ " 691El2; 



With that, we have the three Igusa functions 

h = 2 . 35X12 ^ J2 = 2 -3 3 3^4XJ2 ^ J3 = 2-5. 3^12 +2" 3 -3 2 ^i 2 - (1.2). 
X10 X10 X10 X10 

Igusa shows the equivalence with the definition of these functions in terms of theta 
functions in [10, p. 848]. The analogue of the denominator A = g^ — llg^ appearing 
in (1.1) is the form x\§. The form A is a classical cusp form of weight 12 and xio 
is a Siegel cusp form of weight 10. 

A mathematically natural question is whether we can use formula (1.2) directly 
to evaluate the Igusa functions, thereby bypassing the 6>-functions. The main focus 
of this paper is to give an explicit algorithm to evaluate the Siegel modular forms 
occuring in (1.2) to some prescribed accuracy. Our result gives a relatively easy way 
to analyze the precision necessary for the computation to succeed, and we give a 
rigorous complexity analysis for our method, something which has not been done 
for other approaches. 

Although the asymptotic convergence of our algorithm is slower than the algo- 
rithm using theta functions, our approach has the advantage that there are fewer 
high precision multiplications required in the evaluation, and thus less precision loss 
and fewer rounding errors occur. Furthermore, we give a detailed analysis of the 
Eisenstein series and cusp forms, including an algorithm for computing them using 
classical modular forms of half-integral weight and explicit bounds on the size of 
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the coefficients in their Fourier expansions. Indeed, one of the main contributions 
of the paper is the detailed analysis of various aspects of the computation of Siegel 
modular forms. Finally, our approach may lend itself to improvement in various 
ways and is a new direction in this area which could produce further progress. 
Any Siegel modular form / admits a Fourier expansion 

/(r) = J2 a ( T ) ex P( 2 ™ ^(^Y)) (1.3) 

T 

where T ranges over certain 2 x 2-matrices with coefficients in |Z. We propose 
to evaluate the functions occuring in (1.2) by truncating the sum in (1.3) to only 
include matrices with trace below some bound. The Eisenstein series are Siegel mod- 
ular forms with a considerable amount of extra structure. We show that computing 
the Fourier coefficients of the Eisenstein series ultimately boils down to computing 
Fourier coefficients of classical modular forms of half-integral weight. One of the 
main results of this paper is the following theorem, proved in Section 4. 

Theorem 1.1. For A,Ce Z> and B G Z with B 2 < AAC , the Fourier coefficients 

of the Siegel Eisenstein series E w for all matrices ( b % ' c J satisfying < a < A, 

< c < C , \b\ < B can be computed in time 0((ABC) 1+s ) for every e > 0. The 
constant in the O -symbol depends on the weight w. 

By examining the size of the Fourier coefficients more closely, we derive the following 
result in Section 6. 

Theorem 1.2. Let r G H 2 be given, and let 5 = 5(t) be the supremum of all 
5' G R such that Im(r) — 5'1 2 is positive semi-definite. Assume that 5(r) > 1. 
Assume Xio(t) is non-zero and choose n G Z such that |xio( T )| > 10 _n holds. 
For a positive integer k, let B G Z>o be such that 

POO 

/ 524093t 15 exp(-27rt5(r))dt<10- fc - max{22 ' 6n} (1.4) 

JB-l 

holds. 

Then the following holds: if we approximate the modular forms E4, Eq, xiOi X12 
using their truncated Fourier expansions consisting of all the matrices of trace at 
most B , then the values ji(r), J2(t), J3( T ) computed via the formulas in (1.2) are 
accurate to precision 10 _fc . 

The condition 5(t) > 1 is mostly for esthetic reasons. The proof of Theorem 1.2, 
given in Section 6, readily gives a method to find B in case 8(r) < 1. We assume in 
Theorem 1.2 that we can bound |xio(t)| from below. This lower bound will allow 
us to bound the precision loss that occurs when we divide by Xio(t)- Using the 
explicit bounds on the Fourier coefficients of xio, proved in Section 5, we give a 
simple method to find a value of n in Section 6. This method works in general and 
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does not depend on the value of S(r). Hence, Theorem 1.2 gives an effective method 
to evaluate the three Igusa functions up to some prescribed precision. 

Just as the elliptic j-function is invariant under SL2(Z), the Igusa functions 
ji, j'2, J3 are invariant under the symplectic group Sp4(Z). Hence, we may translate 
the argument r by a matrix M G Sp4(Z) to obtain an Sp4(Z)-equivalent r' G H 2 . 
The value S(t') can be significantly different from 5(r), see e.g. Example 7.1. Before 
applying Theorem 1.2, we therefore move, using e.g. the method from [23], r to the 
'standard' fundamental domain for Sp4(Z)\H2 described in [8]. 

The outline of the article is as follows. In Section 2 we recall basic facts about 
Siegel modular forms and their Fourier expansions. Section 3 introduces Jacobi 
forms and their relation to Eisenstein series. The approach we follow in this section is 
'classical' and most likely well-known to experts working with Siegel modular forms. 
In Section 4 we go one step further, and relate Jacobi forms to classical modular 
forms of half-integral weight. This gives a very efficient method of computing the 
Fourier coefficients of the 2-dimensional Eisenstein series. The functions %io and 
X12 are Siegel cusp forms, and we explain in Section 5 how to compute the Fourier 
coefficients of these forms. We investigate the convergence of the Fourier expansions 
of i?4, Eq, Xio and X12 in Section 6. This leads to the proof of Theorem 1.2. A final 
Section 7 contains two detailed examples. 

2. Siegel modular forms 

Let H2 = {t E Mat2(C) | r = r T ,Im(r) > 0} be the Siegel upper half plane. 

With J — [ _ 1 2 j, the symplectic group Sp4(R) is defined as Sp4(R) = {M E 

GL/4(R)|M JM T = J}. The group Sp4(R) naturally acts on the Siegel upper half 
plane via 

'a b\ ar + b 



c dJ ct + d 

where dividing by ct + d means multiplying on the right with the multiplicative 
inverse of the 2 x 2-matrix ct + d. The matrix — 12 acts trivially, and it is well 
known that the automorphism group of H 2 equals PSp4(R) = Sp4(R)/{±l 2 }. 

A holomorphic function / : H2 — > C is called a Siegel modular form of weight 
w > if it satisfies 

f({ a c b d )r)=det(cT + drf(r) 

for all r and all matrices in the subgroup Sp4(Z) C Sp4(R). The integer w is called 
the weight of the form /. Whereas we have to demand that / is 'holomorphic at 
infinity' for classical modular forms H — > C, this is not necessary for Siegel modular 
forms. Indeed, the Koecher principle implies that / is bounded on sets of the form 
{r G H 2 I Im(r) > al 2 } for a > 0, see [14]. 

The matrix [ 1 ) is contained in Sp4(Z), and a Siegel modular function / is 

invariant under the transformation r 4 t + 1. In particular, a Siegel modular 
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function / admits a Fourier expansion 

/(r) = ^a(T)exp(27rzTr(Tr)). 

T 

Here, the sum ranges over all symmetric matrices T G Mat(|Z) with integer diag- 
onal entries. The coefficients a(T) are called the Fourier coefficients of /. By the 
Koecher principle, they are zero in case T is negative definite. 

We embed the group GL 2 (Z) in Sp 4 (Z) via M H> (™ (m t } -i)- As M t has 
determinant ±1, we see that a Siegel modular function / is invariant under the 
transformation r h-> MtM t for M G GL2(Z). This invariance is the key ingredient 
in the proof of the following well known lemma. 

Lemma 2.1. The Fourier coefficients a(T) of a Siegel modular form f satisfy 
a(M T TM) = a(T) for every M G GL 2 (Z). 

Proof. Writing r = x + iy with x, y G Mat2(R), the Fourier coefficient a(T) is 
given by 



t (T)= f f(r)e- 2 ™ T ^dx. 



Here, dx means the Euclidean volume of the space of x-coordinates and the integral 
ranges over the 'box' — 1/2 < Xij < 1/2. Using the invariance of / we compute 

a(M T TM) = /"/(MrM T )e- 2 ^ Tr (™ rMT )dx, 

and the lemma follows. □ 

In the remainder of this section we investigate how many different values a(T) 
attains for a fixed value of n — det(T) > and a fixed Siegel modular form /. 

To a matrix T — I , % ' J with a,b,c G Z we associate the binary quadratic 

form fx = aX 2 + bXY + cY 2 of discriminant b 2 — 4ac = —An. An explicit check 

shows that for M = ( a ^ ) the quadratic forms associated to M T TM equals 

f M T T M = MaX + (3Y, 7 X + SY), 

which means that the GL2(Z)-action on H2 is compatible with the GL2(Z)-action 
on quadratic forms of discriminant — An. In fact, the GL.2(Z)-action on quadratic 
forms originally considered by Lagrange is not used that much as it leads to a 
'wrong' kind of equivalence. For quadratic forms, the 'correct' action is the action 
of the subgroup SL 2 (Z) C GL 2 (Z) studied by Gaufi. The difference between these 
two actions is implicit in the following lemma. 
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Lemma 2.2. Fix a Siegel modular form f and n G \7i > q. Suppose that —An is a 
fundamental discriminant and let O be the maximal order of Q(y/—n). Then the 
set {a(T) | det(T) = n} has size at most §(#Pic(£>) + #{a G Pic(C) | 2a = 0}). 

Proof. If —An is fundamental, then any integer binary quadratic form aX 2 + bXY+ 
cY 2 of discriminant — An is primitive. The set of PSL2(Z)-equivalence classes of 
primitive quadratic forms of discriminant — An is in bijection with the class group 
Pic(C) via aX 2 + bXY + cY 2 ^ aZ + ~b+V=^ z by [ 2 , Th. 5.2.8]. 

It remains to investige when a GL 2 (Z)-equivalence class decomposes as 2 disjunct 
SL 2 (Z) -equivalence classes. If a fractional CMdeal o is GL 2 (Z)-equivalent but not 
SL 2 (Z) -equivalent to b, then b equals the inverse a -1 and we have 2a 7^ 0. The 
lemma follows. □ 

For the general case of not necessarily fundamental discriminants, we note that any 
binary quadratic form aX 2 + bXY +cY 2 of discriminant — An determines a primitive 
quadratic form (aX 2 + bXY + cY 2 ) / gcd(a, 6, c) of discriminant — An/ gcd(a, 6, c) 2 . 
Arguing as in the proof of Lemma 2.2, we see that the set {a(T)| det(T) = n G 
^Z>o} has at most 

l - J2 #Pic(0) + #{a G Pic(O) I 2a = 0} 
o 
elements. Here, the sum ranges over all imaginary quadratic orders O that contain 
the order of discriminant — An. 

Corollary 2.3. Let m be the index of the order of discriminant —An in the maximal 
order of the quadratic field Q(^—n) and let ip denote the Euler ip- function. For a 
fixed Siegel modular form f , the set {a(T)\ det(T) = n G |Z>o} then has as most 
2y/nlog(An)(m/ip(m)) 2 elements. 

Proof. The class number for the imaginary order of discriminant D is bounded by 
l-D) 1 / 2 log \D\ by [18, Sec. 2]. The result now follows from the class number formula, 
see e.g. [17, Sec. 1.6]. □ 

3. ElSENSTEIN SERIES 

For w > 0, the space M w of Siegel modular forms of weight w has a natural structure 
of a C-vector space. For even w > 4, the primordial example of a degree w Siegel 
modular form is the Eisenstein series E w defined by 

E w (T) = J2(cr + d)- w . (3.1) 

c.d 

Here, the sum ranges over all inequivalent bottom rows (c d) of elements of Sp4(Z) 
with respect to left-multiplication by SL(2, Z). The restriction w > A comes from 
the fact that the expression in (3.1) does not converge for w = 2. 

The direct product M = U^ =0 M w has a natural structure of a graded C- 
algebra. By restricting the product to even w, we get a graded subalgebra M e . The 
following lemma gives the structure of these two algebras. 
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Lemma 3.1. The Eisenstein series E±, Eq, E w and E\2 are algebraically inde- 
pendent and generate M e . There exists a polynomial P in 4 variables such that M 
is isomorphic to M e [X]/(X 2 — P{E^ Eq, E±o, E12)). The element X corresponds to 
a Siegel modular form of weight 35. 

Proof. The first statement can be found in [11, pp. 194-195]. The second statement 
is proven in [10] with an explicit polynomial P at page 849. □ 

The remainder of this section is devoted to deriving a 'formula' for the Fourier 
coefficient a(T) of the Eisenstein series E w . The approach we follow is intrinsically 
related to the theory of Jacobi forms, see [7] for a good introduction. Let / : H2 — > 

C be a Siegel modular form of weight w. We write r G H2 as r = f T * ^ J . Because 

/ is periodic with respect to T2, it admits a Fourier expansion 

00 
/(r)=]TyUn,£)e 2 ™ 2 

m=0 

where ip m is a function from H x C to C. The functions ip m have the following 
properties: 



o 



vWfrS, ^) = (en + dre^™/(^ Vm (Ti, e), ( a c b d )e SL 2 (Z) 

^ m (r,e + Ar + ^) = e" 2 — ( A2 ^+ 2A£ )^ m (r 1 ,e), (A, M ) G Z 2 



ip m admits a Fourier expansion of the form > > c(n, r)e 7I "*( nT i+ r£ )_ 



Z^ Z^ 

n=0 r(EZ 



r < 4 ri rri 



The first two properties follow from the transformation law of Siegel modular forms 
under the symplectic matrices 



(a b 0^ 

10 

c d 

Vo 1/ 



and 



/l // \ 

A 1 fi 

1 -A 

VO 1 / 



and the third property follows from the Koecher principle. 

A holomorphic function j:HxC-yC satisfying the three properties above for 
some w and m is called a Jacobi form of weight w and index m. Jacobi forms can 
be seen as an 'intermediate' between Siegel modular forms and classical modular 
forms. Indeed, the 'Fourier coefficients' of a Siegel modular form of weight w are 
Jacobi forms of weight w and for a Jacobi form g, the function g(r, 0) is a classical 
modular form of weight w. 

The space of all Jacobi forms of weight w and index m is denoted by J w , m , and 
we have maps 

npr 
m>0 
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where pr denotes the projection onto the first factor. For this article, the key prop- 
erty of Jacobi forms is that we can also construct a map J Wj i — > M w which will 
allow us to identify certain Siegel modular forms with its 'first' Jacobi form. As 
we have J Wj i = for odd w by [7, Th. 2.2], we restrict to even weight w for the 
remainder of this section. 

For m > 0, we define the 'Hecke operator' V m : J W; i — ¥ J w , m as follows. For 

9 £ J w ,i with Fourier expansion Y^ n r c(n, r ) e 27I "H nT i+ r£ ) 5 we p U t 



Vm(g) = J2\ E aW ~ lc { 



n , r \ a|gcd(n,r,m) 



nm r\ { 2iri(nT 1 + r£ ) 

9 ' 



for m > 0. This is the natural generalization of the Hecke operators for classical 
modular forms, see e.g. [21, Prop. VII. 12]. For m = 0, we put 

AW y V w n ^ 

with a n (x) the sum of the nth powers of the divisors of x and B w the wth Bernoulli 
number defined by t/(e t — 1) = Y^=oB n t n /n\. In particular, the function Vo(g) is 
a multiple of the classical Eisenstein series of weight w. It is not hard to show that 
the function _ 

^(9)=J2 V ^9)(ri,e)e 2 ^ m ^ 

m>0 

defines a Siegel modular form of weight w, see [7, Th. 6.2]. 

Lemma 3.2. The map \1/ : J W; i — ¥ M w is injective. 

Proof. This follows directly from the fact that the composition 

*■ - n t pr 

m>0 

is the identity. □ 



\P t— r P r 

Jyo,l ' -Mil) C T | I Jw,m * "w, 1 



We stress that the map \1/ is in general not surjective. The image &(Jk,i) is known 
as the Maafi Spezialschar . However, the Eisenstein series E w e M w do occur at the 
image of a Jacobi form. They are the images of the Jacobi Eisenstein series E^ 
defined by the (rather awkward looking) formula E^(r,z) = 

1 v^ V^ / ,n ,„ ( ( , 9 o,t + b , z cz 2 \ \ 

gcd(c, d) = l 

for iu > 4. Here, a and 6 are integers such that ( a c d J is contained in SL2(Z). 
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Lemma 3.3. We have E w = ^(^-E^ 



Proof. It follows from [7, Th. 6.3] that E w is a multiple of *(£#). Both the 
Siegel Eisenstein series E w and the Jacobi Eisenstein series E^ are normalized with 
constant coefficient 1. The lemma follows. □ 

It is now a straightforward matter to compute the Fourier coefficients of the Siegel 
Eisenstein series. The result is the following theorem. 



Theorem 3.4. Let E w be the Siegel Eisenstein series of weight w, and let T 
b % b i 2 ) e Mat(I 



' G Mat(iZ) be a positive semi- definite matrix with integer entries on the 



diagonal. Write D = b 2 — Aac < and let D be the discriminant ofQ(VD). Then 
the Fourier coefficient a{T) equals 1 for a = b = c = and 

^ ^ d-MD/a*) 

d\ gcd(a,b,c) 

otherwise. Here, Bk is the kth Bernoulli number and a is defined by a(0) = 1 and 

where C is Cohen's function defined by 

C(s - 1, D) = L Do (2 - s)J2^(d) (^) cf-V 2s _ 3 (//d), D = D f. 

d\f V J 

Here, ( denotes the Dedekind (^-function, Lr> is the quadratic Dirichlet L-series, 
H is the Mobius function, (7) is the Kronecker symbol and a n (x) denotes the sum 
of the nth powers of the divisors of x. 

Proof. By [7, Th. 2.1], the term ot(D') equals the Fourier coefficient a(n,r) of the 
Jacobi Eisenstein series E^ with D' = r 2 — An. By Lemma 3.3, we have to apply 
the Hecke operators V m to these coefficients. The theorem follows. □ 

A formula for a(T) is also given in Corollary 2 to [7, Th. 6.3]. In this formula, the 
Bernoulli numbers and the (^-function from Theorem 3.4 are missing. 

We see that Theorem 3.4 gives a much better bound than Lemma 2.2 for the 
cardinality of {a(T) \ det(T) = n E |Z} for the Eisenstein series. Indeed, for fun- 
damental discriminants — An, we have only one Fourier coefficient a(T). In general, 
the number of coefficients is bounded by the number of square divisors of —An 
which in turn is bounded by 0(n £ ) for all e > 0. These bounds hold in general for 
functions in the Spezialschar ^(J Wj i) C M w . Indeed, the Fourier coefficients c(n, r) 
of a function g G J w ,i only depend on the value An — r 2 , cf. [7, Th. 2.2]. 
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Corollary 3.5. Let n^ be the numerator of the kth Bernoulli number B^. Then 
the Fourier coefficient a(T) of the Siegel Eisenstein series E w for the matrix T = 

( 6/2 c ) ^ s contained in the set l/(n TO ri2 10 -2)Z C Q. 

Proof. As we have £(3 — 2w) = —B2 W -2/(2w — 2) all we have to do is examine 
the denominator of the value Ld (2 — w) occuring in Theorem 3.4. This is most 
easily done using p-adic L-series as in [3, Ch. 11]. The corollary follows from [3, Cor. 
11.4.3] except in the following case: the discriminant of Q(V^ 2 — 4ac) equals — p for 
an odd prime p with w — 1 = (p — l)/2 mod (p — 1). If this is the case, we a priori 
find that the denominator of the L- value could be divisible by p(w — 1). However, 
the prime p then satisfies (p — 1) | (2w — 2) and by the Clausen- von Staudt theorem 
[3, Cor. 9.5.15] the prime p also occurs in the denominator of L?2 TO -2- Finally, w — 1 
is a divisor of the denominator of £(3 — 2w). □ 

Corollary 3.6. The Fourier coefficient a(T) of the Siegel Eisenstein series E w 

for the matrix T = L a /2 b/ c 2 ) satisfies \a(T)\ = 0({4ac - b 2 ) w ~ 3 / 2 ) if b 2 - 4ac is 
non-zero. 



Proof. Using the functional equation for Dirichlet L-series, see e.g. [3, Th. 10.2.6], 
we bound Ld (2 — w) = 0(Dq ' ). The inequalities 

^M = V-UV-L=C(n)<oo 

d\x d=l 

give a n (x) = 0(x n ) for n > 1. It follows that the c(D') in Theorem 3.4 is of size 
0(D' W ' ). As Y^d\n d w ~ x /d 2w ~ s is finite for w > 4 and n — >■ oo, the corollary 
follows. □ 

Remark 3.7. It is not hard to make the constant c in the O-symbol explicit. One 

can take 

4w(w - 2)l((w - l) 2 C(2w - 3)((w - 2) 



ir w - 1 ((3-2w)B l 



4. Computing special values of L-series 

The hard part in computing Fourier coefficients of Siegel Eisenstein series is com- 
puting the special values of L-series occuring in Theorem 3.4. If the discriminant of 

the quadratic field Q(^/b 2 — 4ac), corresponding to the matrix I b % c ), is small 

these computations can be efficiently done employing generalized Bernoulli numbers 
as we now explain. 

For n > 1, we let Xn be the quadratic Dirichlet character modulo n and define 
the Xn-Bernoulli numbers Bk(Xn) by the expansion 

Er=i Xn(r)te rt _^B k ( X 



e nt _ I A^ & 

k>0 



ET^^i- ( 41 ) 
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The generalized Bernoulli numbers Bk(Xn) equal the ordinary Bernoulli numbers 
Bk for n = 1 and k > 2. 

Lemma 4.1. For n > 1 and w > 2, we have L n (2 — w) = —B w _i(x n )/(w — !)• 

Proof. See [24, Th. 4.2]. □ 

The values B w -i(xn) can easily be computed using the definition (4.1) for small 
w and n. For evaluating the Igusa functions, we are only interested in the values 
w — 4, 6, 10, 12 and by computing .Bii(xn) we get the other values Bg(x n ), ^s(Xn) 
and £?3(x n ) 'for free'. 

To compute the Fourier coefficients of the Eisenstein series E w for large values 
of D = b 2 — 4ac, we clearly need another method. It is a relatively well-known 
fact that Jacobi forms of even weight and index 1 'correspond to' classical modular 
forms of half-integral weight. Explicitly, for a discriminant D < 0, we define 

«» (D) = c(^) C(, "- 1 ' D) 

as in Theorem 3.4. Here, Dq is the discriminant of the quadratic field Q(\AD) and 
/ satisfies -Do/ 2 = D. We put a w (0) = 1, and a w (D) = if D < is not a 
discriminant. The function H w : H — > C defined by 

oc 

#w(z) = ^ a w (-n)g n (g = exp(2niz)) 

n=0 

is known as Cohen's function. 

Lemma 4.2. Let H w be defined as above. Then H w is a modular form of weight w — 
1/2 for the congruence subgroup r (4). 

Proof. See [4, Th. 3.1], or an alternate proof in [13, Prop. IV. 6]. □ 

Remark. The bound a w (n) = 0(n w ~ 3 / 2 ) from the proof of Corollary 3.6 is in nice 
accordance with the general result that the Fourier coefficients of a modular form 
of weight k are of size (9(n fc_1 ). 

As the C-vector space of modular forms of fixed (half-integral) weight is finite 
dimensional, we can easily compute coefficients of H w given a basis for the vector 
space. It is not hard to show that the function 



e(z) = Y J O n =l + 25>" (q 

n£Z n—l 



e 2niz ) 
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is a modular form of weight 1/2 for To (4). The function 



9{z) = 9\z + 1/2) = (l + 2 ^(-l) n Q n2 ) (q = e 



is therefore a modular form of weight 2. Analogous to the proof of [13, Prop. IV. 4], 
it follows that 9 and 9 generate the C-algebra of all modular forms. The main 
advantage of choosing this basis is that 9 is very lacunary. 

Proposition 4.3. The following equalities hold: 
9 7 + 79 3 9 



H A = 
H 6 = 

Hio = 



-9 11 + 229 7 9 + 119 3 9 2 

32 
-438670 19 + 725876# 15 # + 12824886# n 2 + 8845412 W + 107597# 3 # 4 

22459904 



_ 77683# 23 + 212405# 19 # + 38627902# 15 # 2 + 100820362# n # 3 
12 " ~ 159094784 + 

19313951# 7 # 4 + 42481# 3 # 5 

H 



159094784 

Proof. Using Lemma 4.1, we compute the first few Fourier coefficients of H w for 
w = 4, 6, 10, 12. With the obervation that H w equals an isobaric polynomial in 9 
and 9, we have to solve a system of w/2 equations in w/2 unknowns. The theorem 
follows. □ 

This theorem allows us to compute the first iV coefficients of H w in time 0(A^ 1+0 ( 1 )) 
using fast multiplication techniques. This leads to the theorem stated in the intro- 
duction. An important conclusion is that it is much faster to compute L-values 
simultaneously than to compute them individually. 

Corollary 4.4. For A,C G Z> and B e Z with B 2 < 4AC, the Fourier co- 
efficients of the Siegel Eisenstein series E w for all matrices ( , % ' ) satisfying 



6/2 

< a < A, < c < C, \b\ < B can be computed in time 0((ABC) 1+e ) for every 
e > 0. The constant in the O -symbol depends on the weight w. 

5. Cusp forms 

The techniques explained in Sections 3 and 4 allow us to efficiently compute the 
Fourier coefficients of Siegel Eisenstein series. This suffices for evaluating Igusa 
functions, since these functions are rational expressions in E w for w = 4, 6, 10, 12. 
However, the denominators of the Igusa functions have more structure: they are 
Siegel cusp forms. It is a natural question to ask if we can compute the Fourier 
coefficients of xio directly via Jacobi forms. We explain this method in this Section. 
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Let M^ be the vector space of classical modular forms of integral weight w, 
and let M 1 = \J W>0 M^ be the space of all classical modular forms. It is well 
known that we have M 1 = C[E\,Eq], with E^ the classical Eisenstein series of 
weight w, see [21, Cor. 2 to Th. VII.4]. We define the Siegel operator S : M ->■ M 1 
as follows. For a Siegel modular form / : H2 — > C with Fourier expansion /(r) = 
J2t a (^) exp(27riTr(Tr)) we put 

w)=E«(c:))e— , — c;> 

n>0 V 7 Z 

The Siegel operator is a ring homomorphism M — > M 1 , and it maps Eisenstein 
series to Eisenstein series. In fact, for the Eisenstein series E w , it is the composition 

of the maps 

np r 1 

J W ,m » Jw,0 > M w , 

m>0 

introduced in Section 2. 

A Siegel modular form / is called a cusp form if it satisfies S(f) = 0. Equiva- 
lently, / is a cusp form if and only if the Fourier coefficients a(T) are zero for all 
semi-definite T that are not definite. It follows from well-known identities between 
classical Eisenstein series that 

XlQ = -43867 • 2" 12 • 3" 5 • 5" 2 • 7" 1 • 53 _1 (£ 4 £6 - £10) 

and 

X12 = 131 • 593 • 2" 13 • 3" 7 • 5" 3 • 7" 2 • 337 _1 (3 2 • 7 2 E% + 2 • 5 3 E^ - 691£?i 2 ), 

are cusp forms. The constants in ^10 and X12 should be regarded as 'normalization 
factors'. 

Lemma 5.1. The ideal of cusp forms in M e is generated by xio an d Xi2- The ideal 
of cusp forms in M is generated by xich X12 o,nd a modular form X35 of weight 35 
corresponding to X in Lemma 3.1 

Proof. See [12, Th. 3]. □ 

It is well-known that the cusp forms xio an d X12 are contained in the Maafi 
Spezialschar ^(J/^i), the gest of the proof being [7, Th. 6.3]. A Jacobi form g e J Wj m 
is called a cusp form if its Fourier coefficients c(n, r) are zero for 4nm — r 2 = 0. In 
particular, the map 

npr 
m>0 

maps Siegel cusp forms to Jacobi cusp forms. In weight 10 and 12 we have the 
Jacobi cusp forms 

pio.1 = ^(^6^4,1 - ElE 6A ) and ^12,1 = ^((El) 2 E^ - E 6 E 6A ), 

with El = 1 + 240J2 n>0 a 3(n)q n and E\ = 1 - 504 J] n>0 a 5 (n)g n the classical 
Eisenstein series. The factor 144 should again be regarded as a normalization factor. 
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Lemma 5.2. We have xio — ^(— ¥>io,i/4) and X12 = ^(<£i2,i/12). 

Proof. The cusp forms %io and X12 are contained in the Spezialschar and therefore 
occur as images of Jacobi cusp forms. The spaces of Jacobi cusp forms of weight 10 
and 12 are 1-dimensional by [7, Th. 3.5]. Using Theorem 3.4, we compute 

•(OH 

for a Fourier coefficient of xio and 

c(l,0) = -2 

for the corresponding coefficient of </?io,i- The result for xio follows. The computa- 
tion for X12 yields 

a {Co°i))=l and « 1 >°) = 10 - 

The lemma follows. D 

To compute the Fourier coefficients of (^10,1 and (^12,1 we note that the coefficients 
c ipi0 1 (n, r) and c^ 12 1 (n, r) only depend on the value of An — r 2 > 0. Furthermore, 
the functions 

K 1Q = Y^ <Vio,i ( k )<l k and K i2 = Yl c Vl2il (k)q k 

fc>0 fc>0 

are classical modular cusp forms, for the group r (4), of weight 9| and ll| respec- 
tively by [7, Th. 5.4]. 

Proposition 5.3. Let 9 and 9 be as in Section 4- Then we have 

9 15 9 - W X1 9 2 + 39 7 9 3 - 9 3 9 4 



K 10 
K 12 



4096 
59 19 9- 169 15 9 2 + 189 11 9 3 - 89 7 9 4 + 9 3 9 5 



16384 
Proof. Analagous to the proof of Theorem 4.3. □ 

It should come as no surprise that there are no terms 9 19 and 9 23 occuring in 
Proposition 5.3. Indeed, the forms K±o and K12 are cusp forms and therefore vanish 
at q = whereas 9 does not vanish at q — 0. 

Proposition 5.3 allows us to evaluate the Siegel cusp forms %io and X12 at arbi- 
trary points r G H2. For the proof of Theorem 1.2, we need a bound on the size of 
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the Fourier coefficients of xw an d X12 as well. We need an explicit bound, like the 
bound in Remark 3.7. 

The 'Resnikoff-Saldana conjecture' ([20]) 

|a(T)|=0((detT) w / 2 - 3 / 4+e ) 

for the size of a Fourier coefficient a(T) of a Siegel modular cusp form of weight w 
is known to be false in general. At this moment, the best known result is 



\a(T)\ =0((detT) 



™/2-13/36+e\ 



for every e > 0, see [15]. We will prove in the remainder of this section that the 
Fourier coefficients of gio and (712 satisfy 



\a(T)\ =0((detT) 



w/2-l/2+s\ 



and we will make the constant in the O-symbol explicit. The reason that our 
bound is better than the bound in [15] is that our Siegel modular forms lie in 
the Maafi Spezialschar, and this allows us to give a stronger bound. In fact, we 
will show that if the Lindelof-hypothesis is true, the Fourier coefficients are of 
sizeO((detT)™/ 2 - 3 / 4+£ ). 

First we will bound the Fourier coefficients of Kiq and Kyi explicitly. One ap- 
proach would be to adapt 'Hecke's proof [21, Th. VII. 5] for cusp forms. This tech- 
nique would yield a bound of 0(n 475 ) for (710, where we can make the constant in 
the O-symbol explicit. However, our modular forms have considerably more struc- 
ture and we will use a variant of Waldspurger's formula to obtain a better bound. 

The modular forms K\q and K\i have the property that their Fourier coefficients 
a(n) are zero for n = 1, 2 mod 4. As a consequence, see [7, Sec. 6], both functions 
are Hecke eigenforms. For every Hecke eigenform / of weight w — 1/2, Shimura 
constructs, see e.g. [7, Sec. 5] , an integral weight cuspform g of weight 2w — 2 with 
the property that 

\a(D)\ 2 = [ 9 f j } ^-h ^XD.w- l)|ZT- 3/2 (5-1) 

holds for the Fourier coefficient a(D) of /. In formula (5.1), known as Waldspurger's 
formula, (•, •) denotes the usual Petersson inner product, and L(g, xd, s) is the L- 
series associated to g, twisted by the quadratic Dirichlet character xd- 

Lemma 5.4. Let (7is,fi , 22 be the modular forms associated to K w ,Ki2 respectively 
under Shimura's construction. Then we have 

,t 18,9 ™\ < 75634 and %? 2 > 9 ™\ < 1197339. 

(-ft- 10, ^ 10) {K.i2,K 12 ) 
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Proof. As the space of weight 18 cusp forms is one-dimensional, we have gis = A .Eg. 
Likewise, (722 = Ai?J . By formula (5.1) we have 

{918,918} = L(9i8,XD,fy 8!_| D i8.5 
(K 10 ,K 10 ) \a(D)\i 't^ 1 ' 

for every discriminant D for which the Fourier coefficient a(D) of K\q is nonzero. 
Since we can compute the Fourier coefficients of if 10, it suffices to explicitly evaluate 
the L-series at the center of the critical strip. 
Since gi% is a Hecke eigenform, the formula 

/•OO 

(2tt)-T( S )L(<7 18 , s) = / g l8 (iy)y s dy/y 

Jo 

is valid for all s£C. Analogous to the example in [16], we derive the relation 

L(9is, XD, 9) = ^-(2n/\D\r jr (|) c(D)M2nn/\D\) (5.2) 

n=l ^ ' 

for (?i 8 = ^2 n c{n)n~ s . Here, we write 

r°° si 

8 {x) = y 8 exp(-xy)dy = — exp(-x) (l + x + x 2 /2! + . . . + x 8 /8!) . 
Ji x '' 

The right hand side of (5.2) converges exponentially fast, and since we know the 
Fourier coefficients of gis we easily compute the first bound of the Lemma. 

Since the space of weight 22 cusp forms is also one-dimensional, the bound for 
K\2 follows analogously. □ 

Lemma 5.5. For every e > 0, the twisted L-series associated to the cusp forms 
g 18 and #22 satisfy 

\L(g 18 ,XD,9)\<B(e,9)\D\ - 5+s \L(g 22 , X d, 11)1 < B(e, ll)|£>|°- 5+e 

for all discriminants D < 0. Here, B is defined by 

r,/ x 1 L/ ^9 >/ N9 r(?l+l/2 + £) 



2tt T v ' v ' r(n-l/2-e) 



Proof. Let (7 = Ylm a ( TO )? m be either p 18 or (722, and let 2w be the weight of g. 
With A(s, xd) = (^) S T(s + w - l/2)L(s + w - 1/2, g, xd), the twisted L-series 
for g satisfies the functional equation 

A(s,Xd) = A(l -s,xd) 
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for all sgC. We will bound L(s, g, xd) on a vertical line to the right of the critical 
strip, which by the functional equation gives a bound on a vertical line to the left 
of the critical strip. A variant of the Phragmen-Lindelof theorem will then give the 
result. 

We put P(s) = (t^) S L(s + w -1/2, g,XD) and A{m) = c(m)/m™ -1 / 2 . We have 

L(s + w — 1/2, <7ig, xd) = ^2 m m° •> an< ^ the coefficients A(m) are bounded 
by cro(m) = J2d\m -*- °y Deligne's theorem [6, Th. 8.2]. For any e > and any t G R, 
we bound 



|L(l + £ + «;-l/2,^ 18 , XD )|<^L^<^_^_i =C( i + £ 



We get \P(l + e + it)\ < \ D \ 1+£ ^ 1+e ) 2 = C 1 (e)\D\ 1+e . Using the functional equation, 
we bound 



\P(-e + it)\ 



r(i + r -« + «,- 1/2) s 



r(-e + it + «;-l/2) 

<C 2 (e)(l + |t|) 1+2e Ci(e)|£)| 1+e , 

where the last inequality follows from Stirling's formula. 

By the Phragmen-Lindelof theorem, see e.g. [5, Sec. VI. 4], we can bound 

\P(a + it)\ <C(e)(l + \t\) M ^\D\ 1+£ for alia G [-e,l + e], 

where C(e) = max{Ci(e), C2(e)Ci(e)} is the maximum of the two ^-dependent 
bounds on the vertical lines, and M(a) = 1 + e — a takes the values M{—e) = l + 2e 
and M(l + e) — 0. Taking a = 1/2 and t — 0, we derive 

\L{w,g,XD)\<V^C(e)\D\ l / 2+s , 

which yields the lemma. □ 

Corollary 5.6. Let the notation be as in Lemma 5.5. Then, for every e > 0, the 
coefficients cio,i(n) and Ci2,i(n) of K\q and K 12 satisfy 

cio,i(n) < 320^5(5, 9)n 4 - 5+1/2e and c 12 ,i < 3843 v / 5(e, ll)n 5 - 5+1/2e 

/or all n> 1. 

Proof. Substitute Lemmas 5.4 and 5.5 into Waldspurger's formula (5.1). □ 

Remark 5.7. The only room for improvement in Lemma 5.5, and hence in Corol- 
lary 5.6, is in the use of the Phragmen-Lindelof theorem. This theorem yields a 
factor n 1 / 4 in the bound. We can use stronger results to lower the exponent, but it 
is harder to make the constants explicit. If the Lindelof hypothesis is true, then the 
factor n 1 / 4 can be replaced by n £ (but the constant could in theory be not explicitly 
computable) . 
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Theorem 5.8. Let the function B be as in Lemma 5.5, and define B 2 (x) = 
exp(2 1 / x /(x log 2)). Then, for every e > and any r\ > 0, the Fourier coefficients 
aio(T) anda 12 (T) of xio and X12 satisfy 

\a 10 (T)\ < 320B 2 (7])^/B(e,9)(4detT) 4 - 5+1 / 2£+r > 

\a 12 (T)\ < 3843B 2 (77) v /B(£,ll)(4detT) 5 - 5 + 1 / 2e+r >. 

Proof. The Fourier coefficient of xio f° r the matrix T is bounded by 

v^ rq /4detT\ nnn r— v^ (4detT) 4 - 5+£ / 2 

d|(4detT) V " / d|(4detT) 

The sum on the right hand side is bounded by 

(4detT) 4 - 5+£ / 2 J2 1 <5 2 (? ? )(4detT) 4 - 5+e / 2+r ', 

d|(4detT) 

see e.g. [9, Sec. 18.1.]. The proof for ^12 is similar. D 

Remark 5.9. If the Lindelof -hypothesis is true, then we get a bound \a(T)\ = 
0(n w / 2_3 / 4+£ ) for the Fourier coefficients of xio and xi2- This bound is optimal 
in the sense of the Resnikoff-Saldana conjecture [20]. 

6. Speed of convergence 

In section we carefully analyse the speed of convergence of the Siegel Eisenstein 
series occuring in (1.2), and this will yield Theorem 1.2 without too much effort. To 
analyse the convergence of a Siegel modular function we a priori have to consider 
three variables. We begin by showing that it suffices to look at a 'one-dimensional' 
convergence problem. 

The imaginary part Im(r) of a matrix r G H2 is positive definite. Hence, there 
exists 5 G R>o with Im(r) > Sl 2 , meaning that Im(r) — 51 2 is positive semi-definite. 
We define 

6(t) = sup{5' G R I Im(r) > 5'1 2 } 

to be the 'largest' of all these values. With this notation, we have the following 
lemma. 

Lemma 6.1. Let T — f 6 % ' c J G Mat2(|Z) be positive semi-definite and let 
t G H 2 . Then the inequality 

|exp(27rzTr(Tr))| < exp(-27rTr(T)5(r)) 

holds. 
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Proof. We have an equality | exp(27riTr(Tr))| = exp(— 27rTr(TIm(r))). Since T is 
positive semi-definite, we have Tlm(r) > T5(r). The lemma follows. D 

We have 

oo 
(6.1) E w (t) = ^a(T)exp(27rzTr(Tr)) = J^ J^ a(T) exp(27rzTr(Tr)) 

T t=0 TeS(t) 

where S(t) is the set of all 2 x 2 symmetric matrices of trace t with non-negative 
integer entries on the diagonal and half-integer entries on the off-diagonal. The set 
S(t) clearly has at most 2(t + l) 2 elements for which a(T) is non-zero. 



The technique of 'splitting up' the evaluation of a Siegel modular form as in equa- 
tion (6.1) enables us to find a lower bound for |xio(t)|- The idea is that if we 
have 



Y^ a(T) exp(27rzTr (Tr)) 



res(t) 

t<B 



> 10 



Y a(T) exp(27rzTr (Tr)) 



res(t) 

t>B 



(6.2) 



then the value of |xio( r )| is roughly equal to the left hand side of (6.2). Further- 
more, we can apply the upper bound for the Fourier coefficients of xio given by 
Theorem 5.8 to bound the right hand side of (6.2). Taking B = 2 yields the following 
lemma. 



Lemma 6.2. Let 

Mi = 



1 
1 



M 2 



§1 



M, 



1 



and for e, r\ > 0, put M(e, rj) = 320B 2 (ri) >jB(e, 9), where the notation is as in Theo- 
rem 5.8. If, for any e, r\ > 0, we have \c\ > 10 J 2 °° 2M(e, r ^t 11 + e + 2r i exp(— 2irt8(r))dt 
for 

c = -exp(27rzTr(Mir)) - - exp(27rzTr(M 2 r)) - - exp(27riTr(M 3 r)), 

£i 4r "rr 

then we have |xio( r )| > 9/10|c|. 

Proof. Since xio is a cusp form, there are no matrices T G S(0) U S(l) for which 
the Fourier coefficient a(T) of xio is non-zero. The only matrices T G S(2) for 
which a(T) is nonzero are the matrices M\, M 2 , M%. These matrices have Fourier 
coefficients 1/2, —1/4, —1/4 respectively. Hence, c equals the left hand side of (6.2) 
with B = 2. 

Using Theorem 5.8, we bound the right hand side of (6.2) from above by 



10^2£ 2 max (a(T) exp(-27rTr(T)5(r))} 

t=3 6 *■ ' 
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<20 / M(e,7])t 11+e+211 exp(-2irtS(r))dt, 
where we used the 'AGM-inequality' 4det(T) < Tr(T) 2 . The lemma follows. □ 



Remark 6.3. In Lemma 6.2, we can choose any e and r\. The optimal choice 
depends on the value ofS(r). 

Remark 6.4. If the condition in Lemma 6.2 does not hold for any s,6, we can 
look at the contribution of all matrices of trace 2 and 3. If that majorates the 
contribution coming from all matrices of trace 4 and higher, we have found a lower 
bound on |xio( r )|- 

Proof of Theorem 1.2. The Igusa functions are rational expressions in the Eisen- 
stein series E4, Eq and the cusp forms xio and xi2- The proof consists of 2 parts: 
first we analyse the 'loss of precision' that occurs when applying the formulas (1.2). 
Knowing the precision to which to evaluate the four Siegel modular forms, we then 
carefully analyse the speed of convergence of these series. 
Using Corollary 3.6, we bound 

(5.2) \a(T)\ < 19230 Tr(T) 5 

for a Fourier coefficient of E\ in case det(T) is non-zero. For det(T) = and 
Tr(T) ^ 0, inequality (5.2) holds by Theorem 3.4. We conclude that j-E^r)] is 
bounded by 



f°° 80 144 

1 + / 2 • 19230t 5 (t + l) 2 exp(-2nt5(r))dt < 1 + — -^ + -— ■= + 

Jo °{Tr °{V 7 



76 



5(r) 8 5(r) 7 5(r) 6 ' 

and our assumption S(r) > 1 implies that we may bound this by 302. For Eq(t) we 
get the bound \E 6 (r)\ «1 + 93/5(r) 10 < 94. Using Theorem 5.8 with 77 = 1.37 and 
e = 0.28, we derive the bounds |xio( r )| < 3487 and |xi2 ('?") I < 361893 for the cusp 
forms. 

Using these four upper bounds, it is straightforward to check that if we evaluate 
all four Siegel modular forms up to k + 22 decimal digits, then we know the products 
Xi2(t) 5 , -E , 4(t)xi2(t) 3 and Eq(t)xi2(t) 2 occuring in formula (1.2) up to k decimal 
digits precision. Furthermore, we know by assumption that Xio(t) does not equal 
zero. Let n G Z be the smallest n such that |xio( r )| > 10 _n holds. By dividing 
by Xio( 7 ") 6 7 we lose max{0, 6n} digits precision. Hence, if we evaluate all the Siegel 
modular forms occuring in (1.2) up to / = k + max{22,6n} digits of precision, we 
know the Igusa values ji(r), .72 (t), J3( r ) U P 1° k decimal digits of precision. 

We evaluate the Siegel modular functions £4, Eq,xio,Xi2 using the sum (6.1), 
truncated to only include matrices whose trace is below some bound B. It remains 
to give a value for B such that the function values are accurate up to / decimal 



EVALUATING IGUSA FUNCTIONS 21 

digits. As the speed of convergence of the four series involved is slowest for x 12 , it 
suffices to look at this function. Taking rj = 1.45 and e = 0.1, we have 

Y^ ai 2 (T)exp(27rzTr(Tr)) < / 524093t 15 exp(-27rt5(r))dt 

i c q/- + \ J B — 1 



res(t) 

t>B 



and if the integral is less than 10 l then the contribution coming from the matrices 
of trace larger than B do not alter the first I decimal digits. The theorem follows. 

□ 

7. Examples 

In this section we illustrate the techniques developed in this paper by evaluating 
j'i(t) for two choices of r. 

7.1. Example. We detail the evaluation of the Igusa functions ji,J2,J3 at 

2 + hi 13 + 26i \ TT 
G H 2 
13 + 26z 83 + 14H/ 

to 500 decimal digits of precision. The Igusa functions are rational expressions in 
the Siegel modular forms E±,Eq, xio and X12, cf. Section 1. The idea is to simply 
evaluate these series at r to high enough precision and then apply the formulas (1.2). 
We have the rather low bound S(r) > 0.15 in this case. However, for the purpose 
of evaluating Igusa functions, we may replace r by an Sp4(Z)-equivalent matrix r' . 
It is straightforward to check that the matrix 



r 





/ 1 


-2 


-13\ 


hi % \ 
i 6iJ 


-5 1 



-3 

1 


-18 
5 




^ 00 





V 



(r) 



lies in the fundamental domain for Sp4(Z)\H2 as e.g. described in [8]. We have 
5{t') > 4.3. 

To bound |xio( r/ )l from below, we apply Lemma 6.2. With the notation of this 
lemma, we compute c ~ —1.28 • 10 -28 and the value of the integral is roughly equal 
to 2 • 10~ 15 for (77, e) = (1.5,0.1). We see that Lemma 6.2 does not apply directly. 
However, if we compute the contribution c' coming from all matrices of at most 4, 
then we get c' ~ —1.28 ■ 10 -28 ~ c but we now have 



POO 

20 / 35557t 132 exp(-27rf • 4.3)dt w 1.1 • 10" 



34 



We conclude that |xio(t')I ^ s bounded from below by 1.28 ■ 10~ 28 . 

The lower bound on |xio(t')I yields that we lose 6 • 28 = 168 decimal digits of 
precision in the computation of j'i(r'). However, we also easily bound |xi2(t')I — 
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4.37 • 10 -29 . Hence, we gain 5 • 29 = 145 decimal digits of precision by multiplying 
by Xi2(t')°- The 'net loss' of precision is therefore only 168 — 145 = 23 decimal 
digits of precision. 

Putting everything together, we need to evaluate the Siegel modular forms 
E&, Eq, xio, X12 up to 524 decimal digits precision to know the values of the Igusa 
functions up to 500 decimal digits precision. The integral 

451485i 15 - 2 exp(-8.67rf)d£ 

B-l 

is less than 10 -524 for B = 49 and we hence have to consider all matrices of trace 
up to 49. 

To compute the Fourier coefficients for all matrices ( b % ' c ) of trace at most 49, 

we compute the Fourier coefficients of all matrices satisfying 4ac — b 2 < 2401 = 49 2 , 
with the convention that we only take the matrices of trace at most 49 in the case 
of determinant 0. To compute all the coefficients a(T) for E±, Eq, xio and X12 we 
compute the first 2401 terms of the power series 

00 00 A 

X — -V 2 ~ / v — -\ , . 2^ 4 



e = i + 2j2<i n2 and # = (i + 2^(-i; 



n „n 



Using Proposition 4.3, we compute the first 2401 coefficients of the modular forms 
if 4 and Hq: 

=^H A = 240 + 13440q 3 + 30240g 4 + 138240q 7 + 181440q 8 + 362880Q 11 + 0(q 12 ) 



"jPHq = -504 + 44352q 3 + 166320q 4 + 2128896q 7 + 3792096q 8 + O(^ ] 



The coefficients of the forms Hi are the Fourier coefficients of the Jacobi Eisen- 
stein series Ef . Using Proposition 5.3 we compute the first 2401 coefficients of the 
modular forms K\§ and K\i'. 

^Kio = -l/4 ? 3 + l/2 ? 4 + 4 ? 7 -9 ? 8 -99/4Q n +0(Q 12 ) 
±K 12 = 1 /12q 3 + 5 /6q 4 - 22 /3q 7 -llq 8 +425/ Aq^ + Oiq 12 ). 

The coefficients of K\o and K12 are the Fourier coefficients of the Jacobi cusp forms 
<^io,i and <pi2,i- 

Since the 4 Siegel modular forms we are interested in lie in the Maafi Spezialschar, 
the Fourier coefficient a(T) of one of them only depends on the determinant of T and 
the greatest common divisor of the entries of T. We make an array 'encoding' these 
Fourier coefficients as follows. For every positive integer iV < 2401, we compute its 
square free part iVo and write iV = Nof 2 . For every divisor d \ f, we compute and 

store the Fourier coefficient belonging to a matrix T = I , % ' J with Aac — b 2 = N 

and gcd(a, 6, c) = d. For E4 and iV = 16 we get 

[997920, 1239840, 0, 1239840] 
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for instance. For N = we make a list of all positive integers d < X and store the 
coefficients for the determinant zero matrices with trace d. 

The computations so far were independent of the choice of r = ( T £ * J G H2. We 

let gi = exp(27riri), q^ = exp(27rzr2) and ^3 = exp(2iviz) be the 'Fourier variables' 
of the entries of r. We compute and store the values q± = 1, q±, q\ , . . . , q* and 
likewise for q^. For £ we need to compute both the first X powers of Q and Q~ l 
because the off-diagonal entries of the matrices can be negative. 

The precision needed for this computation is easily computed. Indeed, the max- 
imum bound for a Fourier coefficient is roughly 10 21 and occurs for ^12 and a 
trace 49 matrix. As we need to recognize the values a{T) exp(27riTr (Tr)) up to 524 
decimal digits precision, we need to compute q±, qi and qs with 524 + 36 = 560 
decimal digits precision. 

After making these 4 lists, we now simply loop over a = 0, . . . , X , c = 0, . . . , X 
and 6 = 0,..., [v^acj and for the triples (a, 6, c) with b 2 — 4ac < X we compute 
gcd(a, 6, c) and look up the Fourier coefficient in the stored array. 

We implemented this algorithm in the computer algebra package Magma. We 
did not attempt to be as efficient as possible in our implementation. On our 64-bit, 
2.1 Ghz computer it took roughly 1 second to compute ji(r), 32(1"), 3z(t) up to 500 
decimal digits precision. We have 

j'i(r) = 17399743914575167430246482183.29799... 

for instance. The computation of the Fourier coefficients of the Eisenstein series is 
negligible: the bottleneck is the 'loop' over all matrices f b % ' c J satisfying < 

a < X, < c < X, \b\ < [V4a~c\ . 

7.2. CM-example. The evaluation of Igusa functions is a main ingredient in the 
computation of Igusa class polynomials, which is in turn used to construct e.g. 
hyperelliptic curves with cryptographic properties. We illustrate our algorithm by 
recomputing ji(r) for a small CM-point r. 

Let K = Q(v — 5 + v^) be a quartic CM field. The extension K/Q is cyclic 
and K has class number two. Using [23, Algorithm 1], see also [25, Thm. 3.1], we 
compute that 

, _ /2.4060038200i 0.4595058410z\ 

"" V0.4595058410z 1.9464979789 J 

is an approximation to the matrix r representing the abelian surface C 2 /$((9^), 
where $ is a CM-type for K. We will work with a 50 digit approximation to r. 

As shown in [23], the values ji(r) are in fact integers. Hence, we only need one 
digit past the decimal place to recognize them and we take k = 1 in Theorem 1.2. 
The matrix r already lies in the fundamental domain for Sp4(Z)\H2, and we have 
$(i~) > 1.66. Just as in the previous example, Lemma 6.2 does not apply directly. 
Using Remark 6.4, we compute c ~ —5.3 • 10 -12 , where we include all matrices of 
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trace up to 6. The corresponding integral is roughly equal to 1.2-10 -16 for e = 0.1 
and r\ = 1.45. We conclude that we may take n = 12 in Theorem 1.2. 

Just as in Example 7.1, we bound |xi2(t) | < 3.1 ■ 10~ 12 . We lose at most 1 + 6 • 
12 — 5-12 = 13 digits of precision, and we need to know the evaluations of the four 
Siegel modular forms up to precision 10 -14 . The integral 

524093t 15 exp(-3.32yrt)dt 

B-l 

is less than 10 -14 for B — 9. To get all matrices of trace at most 9, we take all 
matrices ( b % h ' 2 ) satisfying Aac — b 2 < 9 2 = 81. We compute 

j'i(r) = 6202728393749.9999 . . . 

which is accurate enough to derive ji(r) = 6202728393750. 

In this example, it turns out that we only needed to look at the matrices with 
Aac — b 2 < 6. The fact that our bound of 81 was much higher can be explained 
as follows. Firstly, our analysis for the precision loss is for a worst case scenario 
and we actually do not lose 14 digits of precision in this example. Secondly, we use 
the same bound for all the Fourier coefficients of the matrices of a given trace t, 
whereas these coefficients actually vary quite a lot. 
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